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Abstract 

We consider the Higgs boson decay processes and its production, and provide a pa- 
rameterisation tailored for testing models of new physics beyond the Standard Model. 
We also compare our formalism to other existing parameterisations based on scaling 
factors in front of the couplings and to effective Lagrangian approaches. Different for- 
malisms allow to best address different aspects of the Higgs boson physics. The choice 
of a particular parameterisation depends on a non-obvious balance of quantity and 
quality of the available experimental data, envisaged purpose for the parameterisation 
and degree of model independence, importance of the radiative corrections, scale at 
which new particles appear explicitly in the physical spectrum. At present only simple 
parameterisations with a limited number of fit parameters can be performed, but this 
situation will improve with the forthcoming experimental LHC data. Detailed fits can 
only be performed by the experimental collaborations at present, as the full informa- 
tion on the different decay modes is not completely available in the public domain. 
It is therefore important that different approaches are considered and that the most 
detailed information is made available to allow testing the different aspects of the Higgs 
boson physics and the possible hints beyond the Standard Model. 



1 Higgs coupling parameterisations 



The discovery of a new resonance at a mass of 125 GeV, announced by both ATLAS and 
CMS in July this year, has opened a new era in particle physics. In fact, the new particle 
has similar behaviour as the one expected for a Standard Model (SM) Higgs boson, and it 
has been observed in different channels, notably the decay into a pair of photons and into 
a pair of massive gauge bosons [ZZ and W + W~ , where one of the two vectors is virtual). 
If confirmed, the discovery of the hard-sought SM Higgs boson would complete the picture 
for the Standard Model. However, intriguing discrepancies have also been observed, like for 
instance an excess in the di-photon rate and the non-observation of di-tau signal events by 
CMS. Such discrepancies are not statistically significant, thus they may disappear as mere 
statistical fluctuations. Even in this case, knowing the couplings of the new resonance is a 
crucial test for the SM hypothesis. Furthermore, many models of new physics, especially 
the ones addressing the problem of the hierarchy in the electroweak symmetry breaking sec- 
tor, predict sizeable deviations in the Higgs couplings. Precise measurements of the Higgs 
properties, therefore, can give precious information on the kind of new physics that Nature 
chose. The information that can be extracted at the LHC is rather limited to a few channels, 
nevertheless it is important to make the best out of it. In the light of this consideration, it is 
important to choose the most relevant parameterisations for the deviations from the SM cou- 
plings, and use them to present the LHC measurements. Higgs coupling parameterisations, 
meant to perform fits on the available data, may follow different approaches which are not 
completely independent. In the following we will discuss a few of them without aiming at an 
exhaustive description. We shall then discuss our suggestion which is particularly motivated 
by testing models Beyond the Standard Model (BSM). 

From the experimental point of view it makes sense to just parameterise Higgs physics in 
terms of observed quantities such as branching ratios and cross-sections. This is for example 
the case of the parameterisation proposed in Ref. jl], where the relevant cross-sections and 
partial decay widths are multiplied by a suitable factor. The advantage of such an approach 
is its simple link to the experimentally measured quantities. On the other hand, with such 
a choice, correlations among the different parameters are not explicit, in particular between 
tree level and loop induced observables. For example, a modification of the couplings to 
W bosons and top, while modifying tree-level branching ratios and cross-sections for the 
Higgs boson, can also affect the loop-level couplings for the Higgs production via the gluon 
channel or the Higgs decay into two photons. This is in principle not a limitation, but in 
practice if one wants to take these correlations into account a different choice of parameters 
needs to be considered. Another point is: what if, instead of a generic fit, one aims at 
discussing limits for a particular BSM model? In that case explicit correlations among the 
fit parameters should be calculated in order to correctly compute the number of independent 
degrees of freedom for the fit. Therefore, a parameterisation which can easily be connected 
to any model of new physics is also useful. Along these lines, we propose an extension 
of the parameterisation in Ref. [2j, where the contribution of loops of New Physics to the 
H gg and H — > 77 modes is explicitly disentangled from the modification of tree level 
couplings, thus removing correlations among the various parameters. Furthermore, the loop 
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contributions are normalised to the top ones, thus simplifying the interpretation in terms of 
new models. 

Another point of view consists in parameterising physics in terms of effective operators. 
Typically this can be done in a consistent way by organising the operator expansion in terms 
of a given hierarchy by including all possible terms up to a given order. An example is given 
by chiral electroweak Lagrangians |3]|4] . This approach has the advantage of inheriting all the 
standard know-how in effective theories, including the calculation of radiative corrections. 
The possible remarks to such an approach is the large number of effective couplings when 
going beyond the lowest order set and the treatment of possible light degrees of freedom 
beyond the SM particles (this possibility is not completely ruled out as particles in the 
same mass range as the SM Higgs boson might still be possible, see for example [5]). A quite 
detailed and clearly written overview of effective Lagrangians for Higgs physics is given in [6] , 
including the treatment of radiative corrections. 

The question of the choice of a parameterisation does not only depend on the particular 
use or preference for a given formalism, but also depends crucially on the number of fit 
parameters with respect to the number of physically independent data channels. In the 
frequentist approach one computes a \ 2 P er degree of freedom within a given choice of 
model (with a given number of independent fit parameters). The choice of a particular fit 
"model" is a matter of a reasonable rule of thumb, as too many fit parameters compared to 
too few experimental data channels give a poor quality fit [j] 

In the following we wish to discuss how the parameterisation of the Higgs couplings to 
gluons and photons proposed in Ref. [2] can be extended to include tree level couplings mod- 
ifications, and how it compares to other parameterisations, in particular the one proposed in 
Ref. [l]. We will try to show that the extension of [2] is especially useful when the results of a 
fit on the parameters are to be interpreted in terms of specific models of new physics beyond 
the standard model. For the tree level couplings, we follow the same parameterisation as in 
Table 2 of Ref. [I], i.e. we introduce a scaling factor in front of the coupling, Kx where X is 
any massive particle of the SM the Higgs couples directly to. The same scaling factor will 
appear in front of some cross sections and partial decay widths. For example: 

°wh — K w a Wh , o Z h — K z a Zh , a tth — K t a m . [i) 

For the partial decay widths: 



WW 



2 r SM F _ 2 r SM F „2t-.SM r _ r 2 r SM / \ 

K W L WW > L ZZ — K Z L zz , 1 bb — K b L bb , 1 T + T ft r l T+T - , . . . [Z) 



For the Vector Boson Fusion (VBF) cross sections, it is imperative to distinguish the two 
production channels with W or Z fusion: 

&VBF — K w a WF + K Z U ZF . [6) 

So far, the parameterisation is the same as in Ref. [l]: the crucial differences arise in the 
treatment of loop induced couplings, as explained in the following Section. 

1 "With four parameters I can fit an elephant, and with five I can make him wiggle his trunk". Rule of 
thumb attributed to John von Neumann by Enrico Fermi, as quoted by Freeman Dyson in "A meeting with 
Enrico Fermi" in Nature 427 (January 2004) p. 297. 
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2 Loop induced couplings: 



Kjgg Cind K/ryry VS. 



Kjg Cind K/ry 



The parameters introduced so far describe tree level couplings of the Higgs. Typically, 
sizeable modification to such couplings are generated by tree level effects from New Physics, 
like for example mixing of the SM particles with heavier states. Modifications of the loop 
induced couplings, however, deserve a different treatment, because they are directly sensitive 
to any new state that may enter the loop. In Ref. [I], a scaling parameter was also introduced 
to describe the new physics effects in the couplings to gluons and photons, namely: 



OggH — K g a ggH > 



99 



Kj 9 l 99 ' 



77 



77 



(4) 



However, both n g and k 7 depend non trivially on the tree level couplings, in particular Ky/ 
and n t , because a modification of the couplings to W and tops would affect the SM loop 
contribution to the couplings to gluons and photons. So, in general, there is a correlation 
intrinsic in this scaling parameter approach. 

On the contrary, in Ref. j2], we proposed an alternative parameterisation of the couplings 
to gluons and photons that can deal with loop corrections from new physics in an independent 
way from the tree level corrections to the couplings to massive SM states. In this way, the 
parameters are not correlated to each other. Note that a parameterisation equivalent to the 
one we present here has been recently and independently used in [7] for studying the effect 
of vector-like fermions. Furthermore, our parameterisation allows to obtain bounds that are 
more easily interpreted in terms of new physics models. The new parameters, K gg and k 77 , 
enter at the level of the amplitude of the loop corrections. In terms of the partial decay 
widths, we have: 



77 



r 



99 



G F a 2 m 3 H 
128v / 2tt 3 

G F a^m 3 H 



k w A w (t w ) + C] 3 



C 9 t -A t {r t ) [Kt + K gg ] + . 



A(r t ) [n t + k 77 ] + . . 



(5) 
(6) 



where the dots stand for the negligible contribution of the light quarks. The coefficients C*7 
and Cf contain the NLO QCD corrections to the SM amplitudes. Aw and A t are the well 
known W and top amplitudes: 



Mr) 
A w (t) 



T 



(r + (r-l)/(r)) 



--(2r 2 + 3r + 3(2r-l)/(r)) , 



T 



(7) 
(8) 



where r 



4^ aIld 




T < 1 

r > 1 



(9) 
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The amplitudes Aw and A t have the property that they rapidly asymptotise to a constant 
value for large masses of the states inside the loop, i.e. for small r. For a Higgs mass of 125 
GeV, we find: 

A w (r w ) = -8.32 , A t (r t ) = 1.37 ; (10) 

where the top amplitude is very close to its asymptotic value A t (0) = 4/3 ~ 1.33. The 
W amplitude however significantly deviates from the asymptotic value A^(0) = —7. The 
QCD corrections to the top loops can also be computed in the asymptotic limit of large 
top mass, which gives a good approximation of the much more complicated mass dependent 
corrections. One finds [8]: 

c? = i--, C* = l + §^. (11) 

71 Z 71 

In the case of the gluon loop, the corrections, which also includes real emission of an addi- 
tional gluon and splitting into a pair of light quarks, reduces to this simple form if a s entering 
the LO amplitude is evaluated at the renormalisation scale fx = e -7 / 4 m# ~ 22 GeV: in other 
words, part of the correction is encoded in the running of the coupling constant. A more 
detailed discussion of the QCD corrections that can be included in this parameterisation will 
be discussed in a following Section. 

The original simplified parameterisation of [2] can be recovered setting the NLO coeffi- 
cients C*7 and Cf to one and k,w — 1, K t — 1. Note that the contribution of n t and Ky/ 
was effectively included into the loop parameters, and we will discuss in more detail later 
in which case this procedure is allowed. We have normalised the contribution of the new 
physics loops to the contribution of the top loop alone. This is important when we want to 
interpret the fit of the parameters in terms of the properties and nature of the new physics 
running into the loop. 

Neglecting the contribution of light fermions, we can draw a relation between the param- 
eterisation [I] and ours: 



i4(Kt,K gg ) = (K t + K gg ) 2 , (12) 

(k W A W {tw) + Cti A t(Tt) [Kt + Kr, 

{A w {T W ) + C]\A t {r t )y 



These formulas show clearly the correlation between the parameters. This correlation is 
absent in our proposed parameterisation. 



3 Data analysis 

The data analysed by the two collaborations (ATLAS and CMS) is presented in terms of 
measured cross sections in the relevant decay channels and for various selection rules: in the 
following, we will make use of as much information as it is available. For instance, for the 
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H — > 77 channel, the measured signal is available for several selection cuts, while for H — >■ 
ZZ — > 2 x and H — >■ WW, only the total measured cross section is available. In the 

case of the H^W 7 channel, the searches are actually performed in various selection channels, 
however the detailed efficiency of each signal region for the various production mechanisms 
is not publicly available. As in our parameterisation the contribution of the production 
channels to the total cross section is modified with respect to the SM, a meaningful fit is not 
possible without detailed information on the efficiencies. For this reason, in the following, 
we will only consider the 77 and ZZ — > 2 x (/ + /~) channels in the fit. 

In general, for each selection channel i of the experiments, the data being fitted is rep- 
resented by the best fit values of the signal strength /2i (defined as the observed number 
of events divided by the expected number of events for a SM Higgs boson) as well as its 
uncertainty <7j. As described in j9], to compare these values to theoretical expectations, the 
signal strength in one channel must be compared to the one calculated in each model, given 
for a specific channel by: 

* (n <)S« £ p <r™e» BR*"' 1 } 

where n l s is the predicted number of signal events in channel i in the studied model, and 
(n l s ) SM that same number in the SM. For each production mode p (theoretical calculation 
at NLO) the efficiency of selection of a channel % (experimental observation) is given by e l p , 
considered to stay the same with new physics. This assumption is in agreement with our 
parameterisation that only includes corrections that do not change the kinematics of the 
events. At last, BRi and BRf M are the branching ratios of the Higgs boson into the decay 
channel corresponding to the selection channel i, for both SM and studied model. 

In the following we will focus on the decay channels H — » 77 and H — > ZZ* — > £ + £~£ + £~ 
in both the ATLAS and CMS experiments, were we gathered the relevant data from several 
documents. We also had to make some assumptions when not all the information was 



available. For the H — >■ 77 channel in CMS, information was given in 10 : in Table 2 of 



the reference, one can find the product a p = (&p M £ p ) / (J2 P ' a p' Me P ') f° r eacn selection and 
production channel, whereas the best fit values as well as the uncertainties can be found 



in another table on the corresponding TWiki 11 , for a Higgs mass of m# = 125 GeV. 
The results from the ATLAS collaboration, on the other hand, can be found in note [12] , 
where the selection channel efficiencies are given in Table 6, and the best fit values had to 
be extracted from Figures 14a and 14b for 7 TeV and 8 TeV respectively. 

For the H — > ZZ — > £ + £~£ + £~ decay channel however, no tables were given by either 
experiment and we had to assume identical cut efficiencies for all production channels. This 
seems to be a reasonable assumption for an inclusive channel like ZZ, however this may not 
be the case if discriminants based on the kinematic properties of the leptons are used, as it is 
the case for CMS, where different efficiencies for different production channels may arise. The 
assumption of universal efficiency allowed us to replace the efficiency-scaled cross sections a % 
by standard cross sections a p , which were taken from the LHC Higgs Cross-section Working 



Group's Website 13 . Best fit values and uncertainties were extracted from Figure 19 in [14], 



and Figure 16a from the note 15 
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Selection channel 




h 


0"z 




h 




7 TeV 


8 TeV 


Untagged 


61 


3.15 


1.82 


68 


1.46 


1.24 


Untagged 1 


88 


0.66 


0.95 


88 


1.51 


1.03 


Untagged 2 


91 


0.73 


1.15 


92 


0.95 


1.15 


Untagged 3 


91 


1.53 


1.61 


92 


3.78 


1.77 


Dijet Tag 


27 


4.21 


2.04 








Dijet tight 








23 


1.32 


1.57 


Dijet loose 








53 


-0.61 


2.03 



Table 1: CMS results in the H — > 77 channel 10 11 . 



Selection channel 


%9 




0% 




h 




7 TeV 


8 TeV 


Unconverted central low pxt 


92.9 


0.5 


1.4 


92.9 


1.0 


1.2 


Unconverted central high pxt 


66.5 


0.2 


1.9 


72.5 


0.8 


1.7 


Unconverted rest low pxt 


92.8 


2.4 


1.6 


92.5 


0.9 


1.4 


Unconverted rest high p?t 


65.4 


10.3 


3.8 


72.1 


1.8 


1.8 


Converted central low pxt 


92.8 


6.2 


2.6 


92.8 


3.4 


2 


Converted central high pxt 


66.6 


-4.4 


1.6 


72.7 


3.5 


2.7 


Converted rest low pxt 


92.8 


2.7 


2.2 


92.5 


0.4 


1.8 


Converted rest high pxt 


65.3 


-1.7 


3 


70.8 


0.3 


2.1 


Converted transition 


89.4 


0.3 


3.7 


88.8 


5.5 


3.3 


2 jets 


22.5 


2.7 


1.9 


30.4 


2.6 


1.8 



Table 2: ATLAS results in the H — > 77 channel 12 . 



Experiment 




Y,otK°oth\pb} 




a 


CMS 


15.32 


2.19 


0.7 


0.4 


ATLAS 


1.3 


0.6 



Table 3: CMS and ATLAS results in the H -> ZZ -> £+£-£+£- channel 14, 15 
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The data used in our analysis are summarised in Table [T] and [2] for respectively CMS 
and ATLAS in the H — » 77 channel, and Table [3] for H — > ZZ — > £ + £~£ + £~ channel in 
both experiments. When combining with our procedure the 77 ATLAS results, however, we 
met a problem with selection channel Converted central high pxt- Whereas for all the other 
channels the combination of 7 TeV and 8 TeV data gave results in agreement with the ones 



reported in Figure 14c in 12 , it was not the case for this channel. We therefore decided to 
ignore it in 7 TeV as well as 8 TeV data: the following ATLAS plots are combination of all 
remaining channels. 

To fit our parameterisation, we used a x 2 method where each computed signal strength 
is compared to its corresponding data best fit value through 

x 2 = E^^- ( 15 ) 

i 1 

The values of x 2 computed this way are then compared to the exclusion thresholds at 68 
and 95% CL for a m degrees of freedom \ 2 distribution, where m is the number of channels 
i minus the number of independent parameters (i.e. the number of independent parameters 
we include in the fit). 



3.1 Simple two parameter fit 

To start with, one can restrict the fit to the two parameters describing the loop couplings 
only. This can be justified in two cases: 

- new physics only enters via loops while corrections to tree level couplings are small; 

- the most sensitive measurements only involve loop induced couplings, i.e. new physics 
effects on tree level couplings are small. 

In the second case, one can absorb the contribution of n t and Kw into the loop parameters: 

Kgg + K t -l, (16) 

K^ + Kt-l+^p(Kw-l)- (17) 

Note for example that the tih coupling is not measurable at the moment, thus n t can always 
be absorbed in K gg and k 77 . Furthermore, a small contribution to kw ~ 1 can generate 
sizeable effects on k 77 due to the enhancement factor ~ —4.6. In practice, present 
experimental results allow to perform a meaningful fit only with a very restricted set of 
parameters. This relation is also in agreement with the expectation that modifications to 
the tree level couplings, kw and K t , are generated by New Physics effects at tree level. 

In the two parameter case [2], the signal strength can be calculated with all given data, 
as the cross sections for each production channel is proportional to the one in the SM. We 



K 



.'/.'/ 



K 



77 
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can therefore write for the H — » 77 channels: 

yi-\-K, gg ) o gg t gg -\- a VBF t VBF -\- a ZH t ZH -\- a WH t WH nUjj 

^ Tt n SMA 7 n SM i T SM f i i ^-SAf i X ~RRSM 

99 99 VBF VBF ' ZH ZH ' W H W H ^77 

i + [(i + k, /(/ ) 2 - 1] a;, / . 



99 / i i "-77 



1 + [(1 + «: 99 ) 2 - 1] flfl™ V l 6 A w {r w ) + 1 

The factor containing the SM branching BR^ 1 takes into account the change in the Higgs 
total width due to the modification to the gluon couplings (while we neglect the effect of 
the photon channel in the total width). For the H — > ZZ — > £ + £~£ + £~ channel in both 
experiments, the efficiencies being supposed to be the same and the branching ratio into two 
Z being unmodified, the signal strength can be expressed as: 

_ + ^g)" + 1 

1+^^" l+[(l + ^) 2 -l]^ff (l9) 

S '9 

Using the published CMS and ATLAS results, the result of the two parameter fit using 
Kgg and k 77 is given in Figure [l] where only the inclusive 77 data is fitted. As expected, the 
allowed region runs along the line A, that labels points where the inclusive di-photon cross 
section is unchanged. However, the thickness of the allowed region depends on the value of 
Kgg, signalling that the data are already marginally sensitive to the production channel (large 
values of K gg significantly enhance the gluon fusion channel with respect to the others). In 
Figure [2] we also perform a simplified fit adding to the previous data the extra information 
coming from the h — > ZZ channel into four leptons. As in this approximation we assume no 
variation in the decay widths in gauge bosons, adding the new data chops off the region of 
large K gg , where the signal would be too enhanced. 

For reference sample points for the following models are indicated: 

- [♦] fourth generation, where the result is independent on the masses and Yukawa 
couplings; 



[*] Littlest Higgs 16 , where the result scales with the symmetry breaking scale /, set 
here to / = 500 GeV for a model with T-parity (there is also a mild dependence on 
the triplet VEV x, that we set x — 0); 



A] Simplest Little Higgs 17 (see Section 4.2), where the result scales with the W 



mass, also set here to mw> = 500 GeV for a model with T-parity; 



[■] colour octet model 18 , where the result is inversely proportional to the mass 
ms = 750 GeV in the example (and also depends on two couplings set here to Ai = 4, 
A 2 = l); 



5D Universal Extra Dimension model 19 , where only the top and W resonances 



contribute and the result scales with the size of the extra dimension (here we set 
ttikk — 500 GeV close to the experimental bound from electroweak precision tests); 
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-2-1012 -2-1012 

kyy kyy 



Figure 1: k' J7 and K gg at the LHC for a Higgs boson with run = 125 GeV. The two solid lines correspond 
to the SM values of the inclusive 77 channel (A), and the vector boson fusion production channel (B). On 
the left panel, the fit using ATLAS data. On the right, the fit using CMS data. Darker (lighter) blue are 
the 1, 2 a limits. 



- [if] 6D UED model on the Real Projective Plane 20 , with = 600 GeV is set to 
the LHC bound (2l|; 

- [•] the Minimal Composite Higgs [22] (Gauge Higgs unification in warped space) with 
the IR brane at 1/R' = 1 TeV, where only W and top towers contribute significantly 
and the point only depends on the overall scale of the KK masses, as the other param- 
eters are fixed by the W and top masses; 

- [▼] a flat {W' at 2 TeV) and [4k] warped (1/R' at 1 TeV) version of brane Higgs models, 
in both cases the hierarchy in the fermionic spectrum is explained by the localisation, 
and all light fermion towers contribute; notwithstanding the many parameters in the 
fermion sector, the result only depends on the overall scale of the KK masses. 

The numerical values of the parameters, including only k% and both K t and K\y in the loop 
parameters, are given in Table [4] for all the models listed above. The values are computed 
using the results in [2], while the 6D UED model has been computed in [23 ] . Note that the 
fit we present here is based on two independent parameters, while in many of the models we 
show all the parameters depend on a single model parameter, typically the mass scale of the 
new physics. Therefore, one should take any effective parameter fit with the caveat that the 
fit must be redone for a specific model taking into account the actual number of independent 
parameters in the BSM model. Note also the in all cases, except the fourth generation, the 
result scales with the mass of the new particles, therefore a point that falls in the exclusion 
region implies that a stronger bound on the mass of the new states is imposed by the Higgs 
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kyy 



kyy 



Figure 2: n'^ and K gg at the LHC for a Higgs boson with run = 125 GeV including the extra information 
coming from the h — > ZZ channel. On the left panel, the fit using ATLAS data. On the right, the fit using 
CMS data. Darker (lighter) blue are the 1, 2 a limits. 

measurements rather that the exclusion of the model. We should therefore think of the 
model as covering a line of points connecting the benchmark point in the Figures with the 
origin (SM point). 

3.2 Three parameter fit 

As a sample and more general application of our parameterisation, we have also performed 
a three parameter fit, by adding to the two previous variables an explicit dependence on the 
gauge boson couplings. For simplicity, we assume that a custodial symmetry present in the 
BSM models imposes Kz = K\y — K v, so that we can treat deviations from the SM couplings 
of both W and Z with a single parameter. This is usually the case in any reasonable model 
of New Physics. 

In the 3 parameter fit, the expression of the signal strengths in the 77 channel can be 
written as 



/1 , ^2 i , 2 ( SM i , SM i , SM i \ RR 

\ L ' h gg) °gg tgg + K V \ a VBF t VBF ' a ZH t ZH ' a WH t WH) O ti. 




(20) 



rrSMA 1 n SM A , n SM A , n SM A f> pSM 

99 99 VBF VBF ' ZH ZH ' W H W H DA 77 

(1 + K gg ) 2 a i gg + k 2 v {\ -o-* ) / k 77 




2 



1 + (4 - l)(BR^ w+zz ) + [(1 + K 9g f - 1] BRSM V ^A w {t w ) + 1 



For the ZZ channel: 




(l + K gg ) 2 a gg + K 2 v J2ot h °oth 



l + (4 



+ [(1 + S,) 2 " 1] BR s gg M ' 



(21) 



10 



Model 


parameter(s) 


Kw — 1 






k'^{k u k w ) 


4 th generation 







2 


2 


2 


Simplest Little Higgs 


m w > = 500 GeV 


-0.009 


-0.034 


0.067 


0.11 


Littlest Higgs 


/ = 700 GeV 

m w = 500 GeV, x = 


-0.05 


-0.11 


-0.014 


0.23 


colour octet 


m s = 750 GeV 
Ai = 4, A 2 = 1 





0.37 


0.17 


0.17 


5D UED 


m KK = 500 GeV 





0.20 


0.034 


0.034 


6D UED (RP 2 ) 


m KK = 600 GeV 
(R 5 = 1.5 R 4 ) 





1.00 


0.84 


0.84 


composite Higgs 


1/R' = 1 TeV 


-0.04 


-0.04 


-0.03 


0.14 


flat brane Higgs 


mw' = 2 TeV 


-0.005 


-0.45 


-0.47 


-0.45 


warped brane Higgs 


1/R' = 1 TeV 


-0.11 


-0.65 


-1.08 


-0.57 



Table 4: Higgs coupling parameters for various benchmark models: in parenthesis we indicate 
if K t and/or kw are included in the definition of the loop parameters. In the second column, 
the mass parameter the corrections are inversely proportional to (eventual other parameters 
are indicated in parenthesis). Here we only consider corrections to k w generated at tree 
level. 



In Figure [3] we present the results of the 3 parameter fit, by slicing the allowed region for 
Ky = 1. This choice is motivated by the fact that most of the models in Table [4] have very 
small deviations in the coupling of the W. We see that the allowed region is very close to 
the one obtained in the 2 parameter fit in Figure [2j 

4 Model interpretations 

In general, the contribution of the new loops to the parameters can be written as [2]: 

C 1 3 

K 77 = ~C^~ ~^ c > Np Q NP 9hNP £np , (22) 

NP l 

K 99 = (J9 P 2C( r Np) 9hNP e NP , (23) 

NP t 

where Cj^p/ C] and C 9 NP / Cf are the NLO corrections normalised to the SM ones for the top in 
the 77-h and gluon-gluon-h triangle loop vertex, N C) np is the number of colour components 
of the new states running in the loop, Qnp is its electromagnetic charge, C{r NP ) is the 
Casimir of the colour representation of the new state (C = 1/2 for a fundamental), and e^p 
encodes the amplitude of the new physics loop. For masses of the state in the loop larger 
than the Higgs mass threshold, the amplitudes quickly asymptotise to a constant number 
that only depends on the spin of the state: this is already true with a very good precision 
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Figure 3: Three parameter fit at the LHC for a Higgs boson with m H = 125 GeV (both 77 and ZZ 
channels). Here we present a slice of the allowed region for ny = 1. On the left panel, the fit using ATLAS 
data. On the right, the fit using CMS data. Darker (lighter) blue are the 1, 2 a limits. 

for the top amplitude. Therefore, e^p can be simply approximated by a number dependent 
on the spin of the new states: 

€np = 1 forfermions; 

£np = -21/4 for vectors; (24) 
ejyp = 1/4 for scalars . 

The coupling of the new states to the Higgs is expressed in terms of the Higgs VEV depen- 
dence of the new state mass: 

v dm NP (v) 

ghNP = ^ • ( 25 ) 

rriNP ov 

Thus, in order to connect the new parameters n gg and k 77 to new physics models, it 
is enough to know the quantum numbers (charge, colour representation and spin) and the 



Higgs VEV dependent mass of the new particles. Eq.s [22] and [23] also allow to easily draw 
the correlation between the two parameters in specific models. For instance, in models with 
a single new particle: 

^77 _ 3 N c ^pQ 2 NP 

Kgg 8 C{r NP ) 

where we have assumed that the QCD corrections to the amplitudes are the same for the 
top and for new physics. For a top partner (i.e. a particle with the same quantum numbers 
as the top, except the spin): 

^77 ^99 ■ \^) 
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In most of the models listed in Table |4j the result for the relevant parameters depends on 
a single mass scale (with just mild dependence on other parameters of the model). Therefore, 
the prediction cannot in principle be directly compared with any 2 or 3 (or more) parameter 
fit. For a correct interpretation, one should re-scale the significance of the x 2 to take into 
account the presence of fewer parameters. However, the results obtained with an effective 
parameterisation, like the ones advocated here, can be used to give an indication if the data 
favours or disfavours the BSM model under consideration. 



4.1 Higher-order corrections 

It is a well known fact that Higgs physics (in the Standard Model as well as in many 
of its extensions) is largely affected by radiative corrections, in particular because of strong 
interactions. However, the exact calculation of those corrections is a priori model-dependent, 
which makes it impossible to include corrections in an effective Lagrangian without adding 
new parameters. In fact, higher order corrections can generate new operators or kinematic 
structures together with a simple re-scaling of the LO operators. 

The parameterisation we propose, by splitting clearly the effect of tree level modification 
of tree level couplings and New Physics loop effects on the loop induced couplings, allows to 
easily add, at least partially, NLO corrections to the calculation. This is not a completely 
consistent procedure but the largest corrections are included. For instance, we have already 
seen that NLO QCD corrections to the top loop contribution to both gg-h and 77-h vertices 
factorise and give a simple multiplicative factor. We can expect that QCD corrections to 
the New Physics loops have the same structure and can be factorised: contributions of this 



type are already included in Eq.s 22 and 23 via the factors C]^ p and C 9 NP . A precise and 



self-consistent calculation of such effects must be carried out in any given specific model. In 
the BSM model predictions in Table |4j however, we assume that the QCD corrections are 
the same as in the SM top loop for simplicity. This is true for the contribution of coloured 
fermions (which is the most common case) because the corrections are independent on the 
mass of the fermion. 

As already mentioned, this parameterisation cannot include NLO corrections that gen- 
erate Lorentz structures different from the LO SM ones. For instance, electroweak 2-loop 
corrections to the gluon fusion cross sections, which are proportional to the coupling of the 
Higgs to W and Z, cannot be included [24]. Other corrections are loops involving both the 
production and decay process: such effects, however, are expected to be small because we 
are considering processes with the production of a Higgs boson in an s-channel resonance. 

Anyway, even in an effective theory approach based on effective operators the procedure 
to compute radiative corrections implies adding new parameters to deal with the new struc- 
tures. Indeed full consistency of the procedure can be guaranteed but this requires inserting 
counterterms [6] . From the point of view of a fit adding new parameters for describing a few 
per cent modification of the vertices is not necessarily an improvement. For example the two 
loop electroweak correction to h — > 77 is less than 2%, while the effect on gg — > h is of order 
5%. 
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4.2 An example: the Simplest Little Higgs model 

The new parameterisation offers an easy interpretation in terms of new physics models. The 
idea is that the contribution of new physics to the parameters will in general scale in a 
simple way with the mass of the new states: thus, any model will roughly cover a straight 
line originating at the origin (i.e. the SM point). By simply measuring the length of the 
excluded line, one can extract the bound on the new physics mass scale. 

As an example, we present here the case of the Simplest Little Higgs model, described 
in Section 3.3 of Ref. [2]. The model contains both a W partner W and a top partner T. 
Moreover, mixing between the two states generates a modification of the couplings of the 
Higgs to the SM W and top. Therefore, there will be a contribution to Kw and K t from the 
modified tree level couplings, and to K gg and k 77 from the W and T loops. From the W 
sector, the contributions are: 

n w = l-\^, « Ty(W ') = g<l I n gg (W') = 0. (28) 
o rn w , ru m W / 

Here mw> is the mass of the W new particle. From the top sector, we have: 

m? 4 miff , mS m? . , 

Kt = 1 + -\ - --f- , K gg {T) = Kjj(T) = §- . 29 

rrij, 3 77%, 

As Kt is not measurable, one can include its effects in the loop parameters (in this case, for 
illustration, we will leave kw in the fit, so we are intending a 3 parameter fit): 



4 mlr , „ „ „,„ 125 ml, _ 1 



99 3 m;, 77 48 m 2 w , 3m 2 w , 

Thus: 

- the parameters of the fit only depend on mw/mw' and scale like l/m^* (this scaling 
is lost in the K g and k 7 parameterisation); 

- the correlation between the three parameters (in this model) is explicit, n gg = — ^k 77 = 
-A{k w -1). 

This shows explicitly that specific model points (or regions) can be put in a global fit only 
for qualitative and illustrative purposes. In order to exclude a particular model at a given 
confidence level the x 2 per degree of freedom should be calculate in the specific model, as in 
general this quantity is model dependent due to the fact that the supposedly independent 
fit parameters of a global analysis can be correlated in that specific model (in the example 
just given the three parameters reduce to only one independent parameter). 

5 Conclusion 

We have discussed a generalisation of the parameterisation proposed in [2] to include tree- 
level couplings and we showed how it can be used for testing and putting exclusion limits 
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on models of new physics beyond the Standard Model. The most important radiative cor- 
rections, involving QCD NLO corrections, can be easily included in this parameterisation. 
We have compared this formalism to other parameterisations, in particular the one proposed 
in pi. We showed that the two, while sharing the same tree-level structure, are fundamen- 
tally different concerning the treatment the the loop-level couplings. In particular while the 
parameters in [I] are inspired from the experimentally measured quantities, our parameteri- 
sation is tailored to investigate BSM models, keeping track of the specific correlations among 
the parameters. It also allows more easily to interpret mass limits and contributions to the 
loops giving the effective Higgs boson vertices. In fact, the parameters we propose are easily 
calculable in extensions of the Standard Model, and we clearly see in the plots that different 
models cluster in specific directions in the n gg and k 77 parameter space. This property is 
due to our choice to normalise the New Physics loop to the SM top one. We also performed 
a 2 and 3 parameter fit of the CMS and ATLAS results in the H — > 77 and H — > ZZ 
channels, showing that the 2 parameter fit already includes all the necessary information 
and is therefore a good approximation at this stage. More precise measurements of other 
channels will necessary require the inclusion of more effective parameters. We hope that 
this work will trigger interest of adding the proposed parameterisation to the existing ones 
in performing experimental fits on data by the ATLAS and CMS collaborations concerning 
the Higgs boson. 
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